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ABSTRACT 

We investigate how strong gravitational lensing in the concordance ACDM cosmology 
is affected by the stellar mass in galaxies. We extend our previous studies, based on 
ray-tracing through the Millennium Simulation, by including the stellar components 
predicted by galaxy formation models. We find that the inclusion of these components 
greatly enhances the probability for strong lensing compared to a 'dark matter only' 
universe. The identification of the 'lenses' associated with strong-lensing events reveals 
that the stellar mass of galaxies (i) significantly enhances the strong-lensing cross- 
sections of group and cluster halos, and (ii) gives rise to strong lensing in smaller 
halos, which would not produce noticeable effects in the absence of the stars. Even if 
we consider only image splittings > 10 arcsec, the luminous matter can enhance the 
strong-lensing optical depths by up to a factor of 2. 

Key words: gravitational lensing - dark matter - large-scale structure of the Universe 
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1 INTRODUCTION 

The ACDM model, the current standard model of cosmo- 
logical structure formation, is based on a flat universe with 
cold dark matter and a cosmological constant. It has been 
shown to fit a wide range of observations, including the prop- 
erties of galaxies, their clustering, the accelerated expan- 
sion inferred from the apparent luminosity of distant type 
la supernovae, the structure of the high-redshift intergalac- 
tic medium, and temperature fiuctuations in the cosmic mi- 
crowave background ( Spergel et al. [2007 1. 

Further tests and constraints on the parameters of the 
ACDM model are obtained from measurements of gravita- 
tional lensing effects, which were first discovered through 
multiple images of distant quasars ( Walsh et al.||1979 1 and 
highly distorted images of distant galaxies at optical (|Lynds 
and Petrosian|[T986 Soucail et al.|[l987 l and radio (Hewitt 
et al. 1 1988 1 wavelengths. The results of the first optical sur- 
veys for multiply imaged quasar s ([Crampton et al.|1992[ |Yee| 



et al. 1993 



Maoz et al.^ ^1993, ^Surdej et al.|1993 1 were used 



to constrain the value of the cosmological constant (Maoz 
[and Rix|p93l |Kochanek||1993| p96l ). Later quasar strong- 
lensing surveys (Inada et al.J|2008|) have been used to con- 
strain its possible evolution (Oguri et al. 2008') . Many recent 
observations of gravitational-lensing effects around galaxies 



(e.g. Mandelbaum et al.||2006b Simon et al.||2007l, and in 



and around galaxy clusters (e.g. Mandelbaum et al. 2006a 
IClo we et al]|2006| [Natarajan et al.||2007| |Comerford et al.| 
,2006| [Massey et al.||2007b[ " are well explained by the dark- 
matter structures predicted from the ACDM model. Sur- 
veys measuring the weak lensing effects of the matter dis- 



tribution as a whole (e.g. 


Semboloni et al. 


2006 


Hoekstra 


et al.||2006| |Simon et al.||2007 |Massey et al.||2007a| |Ben- 


jamin et al.|2007 


Fu et al.|2008l are particularly promising 



for further constraining the parameters of the ACDM model. 
An open question is whether the observed frequency of gi- 



ant arcs (Luppino et al. 111999 Zaritsky and Gonzalez||2003 



Gladders et al.|2003 1 is consistent with predictions based on 



the ACDM model with parameters favoured by other ob- 



servations (e.g. |Bartelmann et al.||1998| Oguri et al.| 


2003 


Dalai et al.||2004 


Wambsganss et al.||2004| |Li et al.| 


2006 


Meneghetti et al. 


2007 


1. 



The efforts currently underway to improve the measure- 
ment of lensing effects need to be matched by a compara- 
ble improvement in the theoretical predictions. According 
to the ACDM model, most of the matter in our Universe 
is dark. The baryonic matter contributes significantly, how- 
ever, to the inner regions of galaxies and clusters, which 
is where strong lensing is observed. Most theoretical stud- 
ies of strong lensing by galaxies ([ Turner et al.||19 84| [Maoz 
and Rix||r9 93: MoUer and Blain 200 1[ [Huterere t al. 2005 
Oguri|2006||Chae||2007||Moller et al.||2007||Ogurret al.|20(38 



e.g.) use analytic profiles to model both the luminous and 



hilbertSmpa-gar ching . mpg . de 



the dark components of the lenses (see Kochanek||2006 for 
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a review). Many studies of strong lensing by galaxy clus- 
ters use profiles obtained from A'^-body simulations for the 



dark matter, but neglect the luminous matter (e.g. Bartef 
mann et al!][l998l |Wambsganss ct al. 2004l|Li et al.p006 l 



There are, however, also studies of the effect of galaxies on 



giant-arc probabilities in clusters. Meneghetti et al. (20001, 
for example, placed galaxies randomly into simulated dark- 
matter clusters and concluded that, although the galaxies 
change certain aspects of cluster lensing, they do not have 
a strong effect on the formation of giant arcs. |Flores et al.| 
(20001 used analytic profiles for the cluster halo and the 



cluster galaxies and found that the galaxies do only mildly 
increase the cluster's cross-section for giant arcs. |Meneghetti| 
et al.l ([20031 and iDalal et all (120041) studied the effect of a 



large central galaxy in a cluster, and found that even a very 
massive central galaxy does not greatly affect the giant-arc 
cross-section. Wambsganss ct al. ( 2008| observed moderate 
effects of baryons on the frequencies of giant arcs and large 
image splittings by using a simple description of baryon con- 
densation to place galaxies in simulated dark-matter halos. 
Puchwein et al.l (12005) and'Rozo et al.l (12006 1 incorpo- 



rated a treatment of the baryonic component into cluster 
formation simulations and studied its influence on giant-arc 
probabilities. Such simulations currently have difficulty pro- 
ducing a galaxy population which matches observation, so 
their results are not simple to interpret. Until this problem 
is overcome, a hybrid approach that embeds a semi-analytic 
treatment of galaxy formation within an A'^-body simulation 
of dark-matter evolution ( Springel et al.|20M" 2005 1 appears 
the most realistic way to address these issues. 

In earlier work ( jlfilbert et al.|[2007[ [paper I| in the fol- 
lowing), we studied the statistics of strong lensing by shoot- 
ing random rays through a series of lens planes created from 
the Millennium Simulation (Springel et al .|2005 l. This very 
large A''-body simulation of cosmological structure formation 
did not explicitly include gas physics such as radiative cool- 
ing and star formation, and the results presented in |paper"l[ 
did not account for the effects of the stellar components of 
galaxies. 

In this paper, we extend the work of jpaper l| to in- 
clude the gravitational effects of the stars in galaxies, as 
inferred from semi-analytic galaxy-formation models imple- 
mented within the evolving dark-matter distribution of the 



Millennium Simulation (Springel et al. 2005 Croton et al 



2006 De Lucia et al. 2006 De Lucia and Blaizot 20071 



These models, which have been adjusted to be consistent 
with a large number of observations, couple star formation 
in the galaxies directly to the properties of the underlying 
dark matter. They currently provide the most accurate way 
to simulate how the dark matter distribution is populated 
with galaxies. 

The Millennium Simulation has higher resolution and a 
much larger volume than simulations used in previous stud- 
ies of the effects of galaxies on strong lensing. On scales 
above the resolution limit, the simulation provides a more 
realistic matter distribution than analytic models. Further- 
more, the galaxy model we use provides the stellar mass of 
the galaxies, computed from the assembly history of their 
dark-matter halos, as well as their positions with respect to 
the dark matter in the simulation. This is a considerable 
improvement over the simpler recipes used to place galaxies 



et al.|2000[[2003| [Dalai et al.|2004| [Wambsganss et al.|2008 1. 

Thus we hope to obtain more accurate results, particularly 
in those regions where both the luminous and the dark mat- 
ter are important for lensing. This includes the inner few 
arcseconds of most lenses, regions which dominate many 
surveys for strong lensing (Myers et al. 2003| [Oguri et al" 



2006 



20071 



jBolton et aI1|2006|rWillis et al.||2006| [Cabanac et al 



Our paper is organised as follows. In Sec. [2| we sum- 
marise the main aspects of our method for shooting a rep- 
resentative ray sample through the Millennium Simulation, 
and we describe how we incorporate the lensing effects of the 
stellar mass in galaxies. In Sec. [3[ we present results for the 
magnification distribution, for strong-lensing optical depths, 
and for lensing cross-sections as a function of lens halo mass 
and projected distance between lens centre and image. In 
particular, we compare the results obtained for dark matter 
alone to those obtained when the stellar mass of galaxies 
is also included. The paper concludes with a summary and 
outlook in Sec. [4| 



2 SIMULATION METHODS 

Our approach for simulating gravitational lensing closely fol- 
lows that of |paper~T[ The reader is referred to that paper 
for a detailed description. Here, we summarise the main as- 
pects of the method and discuss the extensions needed for 
the present work. 

In order to calculate image distortions resulting from 
the gravitational defiection of light by matter inhomo- 
geneities between the source and the observer, we use a 



Multiple-Lens-Plane algorithm (e.g. 


Blandford and Narayan 


1986||Schneider et al.|1992||Seitz et al.|1994||Jain et al.|2000| 


Pace et al. 


2007 


1. Lens planes are introduced transverse 



into dark matter halos in previous work (e.g. by Meneghetti 



to the line-of-sight, and matter inhomogeneities in the ob- 
server's backward light-cone are projected onto them. Light 
rays are traced back from the observer to their source under 
the assumption that the rays propagate unperturbed be- 
tween lens planes, but are defiected when passing through a 
plane. The ray distortions (more precisely, the distortions of 
infinitesimally thin ray bundles) induced by the lens planes 
are calculated from the projected matter distribution on the 
planes. 



2.1 The dark-matter contribution 

We use the particle data of the of the Millennium Simu- 
lation ( jSpringel et al.||2005| ) to generate the dark-matter 
distribution on the lens planes. The Millennium Simula- 
tion assumes a flat ACDM universe with a matter den- 
sity of JIm = 0.25 in terms of the critical density, a cos- 
mological constant with SIa — 0.75, a Hubble constant 
100/ikms~^Mpc~^ with h — 0.73, a primordial spectral in- 
dex n — 1 and a normalisation parameter ag = 0.9 for the 
linear density power spectrum. The simulation followed 10^" 
particles of mass rrip = 8.6 x IO^/i-^Mq in a cubic region of 
side length L — 500/i~^ Mpc comoving from redshift z = 127 
to z — 0. The effective resolution reached near the centres 
of dark-matter halos is comparable to the comoving force- 
softening length of 5h~^ kpc. During the run, 64 snapshots 
of the simulation were taken and stored on disk. 
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Along the line-of-sight, we place one plane for each 
snapshot at the corresponding distance from the observer, 
resulting in 44 planes for sources at redshift 2^ = 5.7 (the 
highest redshift we consider). For each lens plane, the parti- 
cles within an oblique slice of appropriate thickness through 
the corresponding snapshot are projected onto a hierarchy 
of meshes (with a spacing of 2.5h~^ kpc comoving for the 
finest mesh). The projected matter distribution is smoothed 
by an adaptive scheme to reduce shot noise from individ- 
ual particles while retaining a resolution of about 5h~^ kpc 
comoving in dense regions. Fast-Fourier- Transform (FFT) 
methods ( Frigo and Johnson||2005 I are employed to calcu- 
late the dark-matter lensing potential from the projected 
particle distribution on the planes. The second derivatives 
of the potential, which quantify the dark-matter contribu- 
tion to the distortion of light rays passing through the lens 
planes, are calculated by finite difference and bilinear inter- 
polation. 



2.2 The stellar contribution 

The Millennium Simulation did not explicitly simulate the 
physics of star formation. Rather this was done in post- 
processing by applying several semi-analytic models to halo 
merger trees generated from the stored output in order to 
follow the formation and evolution of the galaxies. In this pa- 
per we use the catalogue made pub licly available by |Lemson| 
and the Virgo Consortium] (|2006[ P1 based on the model by 
De Lucia and Blaizot (20071 to obtain the main properties 



of the galaxies. Besides many other quantities, the catalogue 
provides the positions, stellar disc and bulg^ masses, and 
disc radii. Ifere, we restrict our analysis to galaxies with 
Afsteiiar J? 10^/i~^Mq, which wc project onto the same set 
of lens planes as the simulation particles. As a test, we also 
tried a higher mass limit A/stciiar > 10^°h-^MQ, finding 
rather minor changes. Clearly our adopted limit is quite suf- 
ficient to include all significant lenses, at least according to 
this model. 

For each galaxy, we approximate the projected matter 
distribution of the disc component by an exponential surface 
density profile with comoving scale radius rgdisc and total 
mass Afdisc taken from the galaxy catalogue data[^ 



Sdisc('") 



Mdisc 



exp 



dis 



(1) 



Here, Edisc(7') denotes the comoving surface mass density at 
projected comoving distance r on the lens plane. 



The galaxy model of De Lucia and Blaizot ( 2007 1 pro- 



vides bulge masses Mbuige, but not the bulge radii. We thus 
use an empirical relation to calculate the bulge radii from 



^ littp:/ /www. mpa-garching.mpg.de/Millennium 
^ Here, 'bulge' means the spheroidal stellar component of a 
galaxy, e.g., the stellar bulge of a disc galaxy or the all stars 
in an elliptical galaxy with no disc. 

^ The catalogue does not contain Mdiac ^^rid disc explicitly, but 
provides the total stellar mass Msteiian the stellar bulge mass 
Mbuige, the physical disc size radius -Rdisci a^ucl the redshift 2. 
According to the underlying galaxy-formation model, Mdisc = 
A^stcllar - A/bulge and rgdisc = -Rdisc(l + ^)/3. 



the bulge masses. [Shenet aL| ( |2003| f] studied the size distri- 
bution of 140 000 galaxies in the Sloan Digital Sky Survey 
(SDSS) and found a relation Rc cc M^^^^^^ between the me- 
dian physical effective radius Rc and stellar mass Matoiiar of 



early-type galaxies at redshift z < 0.3. Trujillo et al. (20061 
combined resuhs of SDSS, GEMS, and FIRES to study the 
evolution of galaxy sizes between redshift z = and z ~ 3, 
and found Rc(x{l + 2)-°"^ for the mean radius _Rc of early 
types at fixed stellar mass. For the effective bulge radius 
f'o bulge (measured in comoving units), we combine these re- 
lations into: 



J^e bulge = + 



,r X 0.56 

1%) X 0-54/^-^ kpc. (2) 



The bulge component of each galaxy is then approximated 
by a spherical de- Vaucouleurs profile ( de Vaucouleurs|1948 1 : 



-'bulge 



(r) 



94.5Mbuige 



exp 



c bulge 



-7.67 



bulge 



1/4 



(3) 



Here, Sbuige(r) denotes the projected bulge surface mass 
density at projected comoving radius r. 

For each galaxy, we use analytic expressions to calcu- 
late the ray distortions induced by the mass distributions ^ 
and (|3| ( see, e.g. 



Cardone||2004 1. The disc and bulge con- 



tributions of all galaxies on the lens plane are then summed 
to obtain the stellar contribution to the distortion of light 
rays passing through the plane. 

The profile ([T]) for the disc component closely follows 
the projected mass distribution of a stellar disc seen face- 
on. This, of course, neglects the effects of disc inclination 
with respect to the line-of-sight, which on average increases 
the cross-section for strong lensing by isolated disc galax- 



ies (e.g. [Mailer et al.||1997 


Wang and Turner||1997 Bartel- 


mann and Loeb||1998| |Keeton and Kochanok 1998| |M611er 


and Blain||1998l. Similarly, 


the model (|3| for the bulge ne- 



gleets ellipticity. Furthermore, recent observations (Trujillo 
et al.'2007) indicate a stronger evolution of galaxy size with 



redshift for very massive galaxies, leading to smaller radii for 
spheroid-like galaxies with Mgteiiar > lO'^^ M© at redshifts 
z > 1 than our estimate ([2|. Thus our model for the stellar 
component of galaxies may be inaccurate for strong lensing 
occurring close to the centres of galaxies. For strong lensing 
at distances of several effective radii, however, the details of 
the stellar matter distribution do not influence the lensing 
properties very strongly. 

To test the dependence of our results on the partic- 
ular choice of the mass profile for the stars in galaxies, 
we replaced the de- Vaucouleurs profile ([sjl by a Hernquist 
profile (Hernquist 19901, and by a Plummer profile with 
Sbuigo(?') oc (rebuige + T^)~^ and the same effective ra- 
dius Te bulge. For 'cuspy' de- Vaucouleurs and Hernquist pro- 
files, the optical depths and cross-sections for strong lensing 
by both dark and luminous matter agree within ~ 20%, 
whereas the results for the 'cored' Plummer profile were 
noticeably smaller. Hence, our results are not very sensi- 
tive to the details of stellar mass distribution as long as we 
use cuspy profiles for the spheroid stellar component of the 
galaxies. 

The particles in the Millennium Simulation are treated 



^ See also Shen et al. ( 2007 1 for an Erratum 
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as coUisionless particles by the simulation and are used in 
our calculations to obtain the dark-matter distribution on 
the lens planes, but they represent the total mass in the 
simulated part of the universe. One could argue that the 
stellar mass we add to the lens planes should be removed 
somewhere else to balance the mass budget. Here, we re- 
frain from doing this for two main reasons: (i) The mass in 
stars (1% of the total mass aX z = Q) is small compared 
to the mass in collapsed objects (50% at 2 = 0). (ii) Al- 
though the baryons that produce the stars originate from 
within the halos around the galaxies, gas physics increases 
the dark-matter density in the inner part of the halos com- 



pared to coUisionless simulations (e.g. Barnes and White 
1984[ [Blumenthal et aL]|1986[ [Gnedin et al.||2004| ). OnlyTn 
the halo outskirts, which are not relevant for strong lensing, 
is there a net decrease of the matter density. The effects of 
the baryons on the dark-matter profile and their relevance 
for lensing (see, e.g., Puchwein et al. |2005: Jing ct al. 2006j 



Rozo et al.|[2006l |Wamb"s gauss et al.||2008| ) are not the pri- 
mary focus of this work, so we neglect these effects in the 
following. 



2.3 Sampling image distortions 

The light rays traced back from the observer through the 
lens planes to their source define the lens mapping G 1-^ /3, 
which relates the 'observed' angular position of a ray in the 
image plane to the 'true' angular position /3 of its source 
in source plane at a given redshift z^. The distortion 
matrix A — {df3/d9), i.e. the Jacobian of the lens mapping, 
can be calculated from the distortions that the ray bundles 
experience when passing through the lens planes. 

For sufficiently small sources, the distortion matrix 
quantifies the image distortions induced by the defiections. 
In this case, the image magnification fi is given hy jj, — 
(det A) and the eigenvalues Ai and A2 ofA(w.l.o.g. |Ai| > 
IA2I) determine the major-to-minor axis ratio r = |Ai|/|A2| 
of the elliptical images of circular sources ( Schneider et al. 



19921 



To quantify the frequency of images with a given prop- 
erty p, e.g. a large magnification or large length-to-width 
ratio, we define the optical depths 



and 



d^e \m\-'^viQ) 



(4) 



(5) 



where lp(0) = 1 if the image at position Q has property 
p, and lp(0) — otherwise. The optical depth estimates 
the fraction of images with property p assuming a uniform 
distribution of images in the image plane, as is the case, 
e.g., for magnitude-limited surveys and a source population 
with integral luminosity function inversely proportional to 
the threshold luminosity. The optical depth estimates the 
fraction of images with property p assuming a uniform dis- 
tribution of sources in the source plane, which is appropriate 
for volume-limited surveys. 

In order to estimate these optical depths, the dark- 
matter and the stellar contributions to the distortion are 
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Figure 1. The probability density pdf'(/j) of the magnification 
ii for sources at redshift = 2.1. The case of lensing by dark 
and luminous matter (solid line) is compared to the case of dark 
matter alone (dashed line). 



calculated at 160 million randomly chosen positions on each 
lens plane. For every considered source redshift, the distor- 
tions from all lens planes with smaller redshift are combined 
at random to generate the distortion matrices for rays along 
1.6 X 10* random lines-of-sight. Combining the distortions 
from different lens planes at random saves us from comput- 
ing the actual path of light rays, while it allows us to sample 
the image distortions on a very large area under the assump- 
tion that the mass distribution projected on different lens 
planes is uncorrelated. 

The measured fraction of rays with a certain property, 
e.g. a large magnification, is then used as a Monte-Carlo 
estimate for the corresponding optical depth to the cho- 
sen source redshift. The optical depth is obtained by ad- 
ditionally weighting all rays by their inverse magnification. 
For comparison purposes, we repeat the procedure using the 
dark-matter contribution alone, and using the stellar contri- 
bution alone. 



3 RESULTS 

3.1 The magnification distribution 

By binning the magnifications of our random rays, we esti- 
mated the probability density function 



pdf'(M') 



d I 



dfi' 



(6) 



In Fig. [T] the distribution including the stellar contribu- 
tion is compared with the dark-matter only distribution for 
sources at redshift z^ = 2.1. There is little difference be- 
tween the two distributions at magnifications « 1. Only 
in the high-/i tail (containing the strongly focused rays) and 
for very low /i (containing the overfocused rays) do the two 
distributions differ significantly. The increased pdf at high 
magnifications hints at a higher optical depth for large mag- 
nifications when the effects of the luminous matter are in- 
cluded. Qualitatively the same behaviour is found for all 
considered source redshifts. 
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3.2 Strong-lensing optical depths 

From our ray sample, we calculate 

• the fraction of rays witli det A < (called type II), 

• the fraction with det A > and tr A < (type III), 

• the fraction with det A < or tr A < 0, i.e. the sum of 
the two previous classes (type II V ni|^ 

• the fraction with a length-to- width ratio r > 10 for im- 
ages of sufficiently small circular sources, and 

• the fraction with magnification |^| > 10. 

The corresponding optical depths Tp(z^) and Tp{z^) are plot- 
ted in Fig. [2] as functions of the source redshift z^. The op- 
tical depths that account for the stellar mass in galaxies are 
significantly larger than those which neglect it. The effect 
is particularly large at low source redshifts, where the dark 
matter alone is much less efficient in producing strong lens- 
ing. 

In contrast to lensing by dark matter alone, the optical 
depths for > 10 and r > 10 are quite similar when the 
luminous matter is included, the optical depth for im- 
ages of type II is no longer much smaller than r/i, and the 
optical depth Tj^j is even larger than r/n. This implies that 
most images of type III are strongly demagnified. These are 
the images of multiply imaged sources that often remain un- 
detected in observations because they are demagnified and 
close to the bright centre of the lens galaxy. 

The optical depths for joint lensing by luminous and 
dark matter show similar behaviour to cored isothermal 
spheres. For singular isothermal spheres, one of the two 



eigenvalues of the distortion matrix A is unity (Schneider 
et al.|1992 1. Hence, — r, and the cross-sections for r > 10 
and > 10 are equal. Singular isothermal spheres do not 
produce images of type III. If the central singularity is re- 
placed by a small core, the lens acquires a finite cross-section 
for strongly demagnified images of type III, and the cross- 
sections for r > 10 and > 10 separate slightly. Thus, the 
addition of baryons apparently makes the mass profiles of 
strong lenses resemble isothermal spheres with small cores. 

The optical depths t^>io for rays with magnification 
/X > 10 (shown in Fig. [2|3) are very sim ilar to the optical 
depths published by Wambsganss et al. (20081. The optical 



depths for lensing by dark matter alone agree within ten 
percent. However, our results show a larger enhancement 
due to baryons (30-60 percent) compared to their work (20- 
30 percent). The larger values we find may be due to the 
higher effective resolution of our method. 



3.3 Lens properties 

As discussed in |paper"l] for lensing by dark matter alone, the 
properties of most strongly lensed rays, i.e. rays with det A < 
0, tr A < 0, > 10, or r > 10, are predominantly caused by 
a single matter clump along the line of sight. This is equally 
the case if galactic baryons are included as lensing material. 
In order to find these matter clumps, which we refer to as the 
lenses of the rays, we use the method described injpapeFT] 
We determine for each strongly lensed ray the lens plane that 
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Figure 2. Comparison between the optical depths caused by 
dark and luminous matter (solid lines) and the optical depths 
caused by dark matter alone (dashed lines), either assuming a 
uniform distribution of images in the image plane (a) , or a uniform 
distribution of sources in the source plane (b) . Shown are optical 
depths for images of small circular sources of type II (black lines), 
of type III (green), with magnification > 10 (blue) and with 
length-to- width ratio r > 10 (red). 



is sufficient to produce the relevant property in the single- 
plane approximation. Depending on source redshift and the 
property considered, this simple criterion identifies exactly 
one lens plane for at least 60% (usually > 80%) of the rays. 
The redshift of this plane is taken as the lens redshift for 
the ray. 

The resulting lens redshift distribution for rays of type 
II V HI is illustrated in Fig. [3] where the cross-section 
dTii\/ii\/ dz^ is plotted as a function of lens redshift z^ for 
various source redshifts 2^. (Not shown are the lens redshift 
distributions for rays with > 10 and with r > 10, but 
these are very similar.) The inclusion of the stellar mass 
slightly increases the typical redshift of lenses (cf. fig. 4 in 



paper II, but most of the lenses still have z < 2.5 even 



for high source redshifts. The lack of lenses at high redshift 
refiects the lower abundance of massive galaxies and halos, 
as well as the less favourable geometry for lensing at these 
redshifts. 

Following the method in [paper l[ we identify for each 
strongly lensed ray the dark-matter halcj^ associated with 
the lens by locating on the sufficient plane the halo with 



® In all situations relevant for this work, images of type II and 
III, and hence all images of type II V III, belong to sources with 
multiple images. 



^ The dark-matter halos considered here were identified in the 
simulation by applying a friend-of-friend group-finding algorithm 
to the dark-matter particle distribution. 
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the highest ratio Ad/b of its virial mass M (defined as 
the mass within a sphere with mean enclosed density 200 
times the cosmological mean) to the projected distance b 
of its centre from the position where the ray intersects the 
plane. The resulting distributions of halo masses for rays 
of type II V IIL with \fj.\ > 10, and with r > 10 are com- 
pared in Fig. ka, where the cross-sections dr^/dlogM^ are 
shown for = 5.7 as a function of lens halo mass M^. 
For > 10^^/i~^ M0, the cross-sections increase by at 
most 40% when the effects of the galaxies are included. For 

< lO^^^/i"^ Mq, the difference is much larger, however. 
In particular, the mass distribution of strongly lensing halos 
extends to significantly lower values when the galaxies are 
included (see Fig. Bp). There is both a maximum at masses 

« 3 X lO^^/i"Hl0 and a low-mass 'tail'. 

For lower source redshifts = 1.1, there is both a low- 
mass 'tail' and a maximum at masses ~ IQ^^/i"^ M0 
for rays with |//| > 10 (see Fig.[5|. For rays with r > 10 and 
rays of II V III, the distribution becomes bimodal with an 
additional maximum at 2 x 10^^h~^ Mq. 

As can be seen in Fig. |4^, drl^iQ/dlog and 
drl^\>io/9iogM^ are very similar for < lO^^/i"^ Mq. 
Strong lenses with these masses clearly show the lensing 
characteristics of isothermal spheres. Lenses with masses 
> lO^^/i"^ M0, however, have much smaller cross- 
sections for r > 10 than for \fj,\ > 10, which can be inter- 
preted as 'convergence-dominated' lensing, with shear being 
much less important than convergence. 

For halos with mass < lO^^/i"^ M0, the dark mat- 
ter alone is unable to produce strong lensing. However, the 
stellar mass of the galaxies changes this. As can be seen in 
Fig.[4]3, the cross-section is maximal for ~ 10^^ Mq 
if only the stellar mass is considered. Even though such 
galaxies typically have stellar masses Metcllar < lO^ft-^M©, 
and Einstein radii < 1 arcsec, their high abundance 
apparently outweighs the small cross-sections of individ- 
ual galaxies. Galaxies with such small Einstein radii pro- 
duce strong lensing in our simulation because we use an- 
alytic expressions for the stellar contribution to the light 
deflection, which are not subject to the (larger) resolution 
limit of the meshes used for the the dark-matter contribu- 
tion. One should, however, keep in mind that Einstein radii 
te < 1 arcsec for galaxies with A/stdiar < 10^^/i~^ M0 im- 
ply that our approach (which essentially assumes point-like 
sources) may overestimate the strong-lensing cross-sections 
for masses < 10^^h~^ Mq and sources with comparable 
angular extent, e.g. distant galaxies. 

The cross-sections for strong lensing by both dark and 
luminous matter shown in Fig. decrease rapidly for 

< lO^^/i"^M0 and vanish for < 5 x IO^^/i'^Mq. 
Apparently, halos with M < 5 x 10^°h~^ Mq do not con- 
tribute to strong lensing even when the baryons in stars are 
taken into account. 



3.4 Images at larger radii 

The inclusion of the stellar mass in galaxies greatly increases 
the optical depths for strong lensing. As discussed in the 
preceding section, a large part of the increase can be traced 
back to galaxies that have small Einstein radii rE < 1 arcsec 
and hence produce only small image splittings. This suggests 
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Figure 3. The cross-section cJtjj^jjj/c'z'^ for rays of type II V III 
as a function of lens redshift for sources at redshift 2:^ = 1.1 
(solid line), = 2.1 (dashed line), and 2^ = 5.7 (dotted line) 
considering lensing by both dark and luminous matter. 
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Figure 4. The cross-section dr^ /dlog as a function of the 
lens halo mass (see text) for sources at redshift = 5.7. 
Panel (a) compares rays of type II V III (solid line), rays with 
\fi\ > 10 (dashed line), and rays with Ir] > 10 (dotted line) for 
lensing by both dark and luminous matter, (b) compares rays 
of type II V III for lensing by both dark and luminous matter 
(solid line), for lensing by dark matter alone (dashed line), and 
for lensing by luminous matter alone (dotted line). 
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Figure 5. The cross-section dr^/dlogM^ as a function of the 
lens halo mass for sources at redshift = 1.1. Compared are 
the rays of type II V III (soUd Une), rays with \fi\ > 10 (dashed 
line), and rays with |r| > 10 (dotted line) for lensing by both dark 
and luminous matter. 




Figure 6. The optical depth ''"jjyjjj(JS 9) for rays of type II V III 
with lens impact parameters 8. Compared are the depths for 
lensing by both dark and luminous matter (solid line), by dark 
matter alone (dashed line), and by luminous matter alone (dotted 
line) for sources at redshifts = 5.7. 



that the effects of the luminous matter on the optical depths 
is smaller for larger image splittings. 

To measure the impact parameter of a strongly lensed 
ray with respect to the lens centre - and thus obtain an 
estimate for the expected image splitting - we now consider 
the dark-matter subhaloQ and - in the case of lensing by 
both luminous and dark matter - the galaxie^as individual 
lens candidates. For each strongly lensed ray, we identify the 
galaxy or dark-matter subhalo on the 'sufficient' plane that 
has the largest ratio M/b^ of its mass M and the square 
of the projected distance between its centre and the ray. 
This ratio compares the impact parameter b to the Einstein 
radius oc -^/M of the lens candidate under the simplifying 
assumption of a point mass|^ 

The resulting optical depths t/jviii(^ 9) for rays of type 
II V III at angular separations from the lens centre are 
plotted in Fig. [g] for sources at redshift = 5.7. The dif- 
ference between the optical depths for lensing by dark and 
luminous matter and lensing by dark matter alone are very 
large for images at radii < 1 arcsec. However, the optical 
depths become more similar for larger radii so that the dif- 
ference is less than 15 per cent for radii > 5 arcsec, and 
becomes negligible for radii > 10 arcsec. The other optical 
depths for strong lensing exhibit a very similar behaviour 
as a function of minimal distance from the lens. These find- 



ings agree well with results of Oguri ( 2006 1, who studied the 



as identified by SUBFIND ( Springel et al.|2001 1 
* with their stars and their associated dark-niatter subhalo (if 
there is one) 

^ Strongly lensed rays at larger impact parameters w.r.t. the 
halo centre are often lensed by substructure, which produce image 
splittings much smaller than the distance from the halo centre. 
The number of strongly lensed rays at large impact parameters 
w.r.t. the halo centre is therefore rather a measure of the amount 
of substructure present in dark-matter halos than a measure for 
large image splittings. If we identify the mass with the largest 
ratio M/b as lens (as in the previous section), too often a cluster 
main halo is selected, although the lensing is due to a smaller 
non-central galaxy in that cluster (as is revealed by inspection of 
individual cases). For this problem does not occur. 



image-separation distribution using analytical mass profiles, 
Meneghetti et al.| ([2000), Mcncghe tti et al.| ( |2003| and 



[Dalai et al.| ( |2004[ ), who considered the influence of galaxies 
on the probabilities for giant arcs and arclets in clusters. 
Furthermore, our results for the optical depths t^>io(SS 9) 



(not shown) agree well with the optical depths of Wambs- 
[ganss et al.] ( [2008[ ). 



The cross-section 9Tj\viii,>e/91og Af"^ for images of 
type 11 V III at radii ^ 9 is plotted in Fig. [7[ as a function 
of the lens halo mass . If all images are considered, the 
cross-sections for lensing by both dark and luminous mat- 
ter differs quite strongly from the cross-section for lensing by 
dark matter alone. Restricted to images at radii 9 > 1 arcsec, 
however, the cross-sections become quite similar. In partic- 
ular, the cross-section for strong lensing at radii ^ > 1 arcsec 
is very smaU for lenses with < lO^^/i"^ Mq even if the 
stellar mass is included. This confirms our assumption that 
the luminous and dark matter in smaller halos produces 
strongly lensed images only at radii < 1 arcsec. For radii 
^ > 5 arcsec, there are virtually no lenses with masses be- 
low AI^ < 2 X IO^^Mq/i^, and the cross-sections are almost 
identical. 

The optical depths r^(z^) for strong lensing by both 
dark and luminous matter as a function of source redshift 
are compared to the optical depths caused by dark mat- 
ter alone in Fig. |8[ for images with lens impact parameters 
5 arcsec. These optical depths can be interpreted as being 
restricted to image splittings > 10 arcsec. The enhancement 
due to the stellar mass is larger for sources at lower redshift. 
For z^ — 1, the stellar mass increases the optical depths by 
up to a factor two. For larger source redshifts the increase 
is < 1.5 depending on source redshift and the property con- 
sidered. Compared to the other optical depths, the optical 
depth for \fi\ > 10 is affected least by the stellar matter. 
Note that the difference in the optical depths for lensing by 
both dark and luminous matter and by dark matter alone 
is larger for than for because of additional differences 
in the magnification distribution. 
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Figure 7. The cross-section Stjivih >e/9logM^ as a function 
of the lens halo mass (see text) for lensing by both dark and 
luminous matter (solid line) and lensing by dark matter (dashed 
hne) for sources at redshift = 5.7. Compared are the cross- 
sections for rays of type II V III regardless of their impact param- 
eters (black lines), at impact parameters 9^1 arcsec (blue lines), 
and at 9 ^ 5 arcsec (red lines). 
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Figure 8. The optical depth as a function of source redshift 
for images of small circular sources of type II V III (black lines), 
with large magnification (blue) and with large length-to-width 
ratio (red) at lens impact parameters 0^5 arcsec. Compared are 
the optical depths for lensing by both dark and luminous matter 
(solid lines) and by dark matter alone (dashed lines). 



3.5 Quasar lensing 

The optical depths for images lensed by luminous and dark 
matter at radii ^ 5 arcsec (Fig.[8| are an order of magnitude 
lower than the corresponding depths for images at radii ^ 
0.5 arcsec (not shown). This ratio is consistent with results 
by Inada et al. (20081, who identified ten lensed quasars 



with image splittings 1-2 arcsec and one lensed quasar with a 
splitting of 15 arcsec in a statistical sample of 22 683 quasars 
with redshifts 0.6 ^ 2 ^ 2.2 from the Sloan Digital Sky 
Survey (SDSS). 



The quasar sample of Inada et al. ( 2008 1 has been used 



by Oguri et al. ( 2008 1 to constrain the matter content and 
dark-energy evolution of the Universe. [Hennawi et al.| ( |2007[ ) 
estimated the number of expected number of quasars with 
very large image separations in similar samples. To obtain an 
estimate of the expected total number of lensed quasars in 



the sample of Inada et al. ( 2008 1 from our data, we assume 



that every strongly lensed source with image splitting 
1 arcsec has one image of type II at radii 0.5 arcsec from 
the lens centre. The quasar sample is limited to apparent 
i-band magnitudes rrii < 19.1. For the calculation of the 
magnification bias, we use the quasar luminosity function 



obtained by Richards et al. ( 2005 1 from the 2dF-SDSS LRG 



and Quasar Survey and Groom et al. (|2004 1 

$*(2 



La 



+ 



La 



L^iz) 



(7) 



Here, Lg denotes the rest-frame (7-band luminosity, z the 
quasar redshift, and ^"{z) a normalisation constant. The 
bright-end slope a = —3.31, and the faint-end slope /3 = 
— 1.45. The break luminosity L*{z) is parametrised as 



Ll{z) 



-|^gfclZ + fe22^-0.4M*{0)^ 



with M*g (0) 



-21.61, fci = 1.39, fca 



so (8) 
0.29, and Lgo as a 



luminosity standard. 

We then convert the apparent i-band magnitude limit 
m''™ = 19.1 to the corresponding g-band luminosity 



(9) 



where DM(2) denotes the distance modulus to redshift 2, 
and the iC'-correction Ki^g{z) 'corrects' between observer- 
frame i-band and rest-frame g-band. We 
correction discussed in [Richards et al. (20061, 



use the K- 
i.e. we as- 



sume Ki^g = Ki^ii{z) — (3.664 with Ki^ii{z) given by table 4 



Richards et al. ( 2006 1 



The quasar lensing cross-section o-ql(2) as a function of 



redshift z is then calculated from our ray sample by ( Schnei- 
der et al.|19"92t : 



0-Ql(2) 



5Z lM(i;2:)r^lQL(i 

i = l 



^{Lg-z)dLg 



/'OC 



,^z)HLs\z)ALg 
z)r^'^Qh{i\z)x 



(10) 



Fi.W;z)\ 



\z)) 



F(x-"™(2)) 

Here, Arrays denotes the total number of rays, z) de- 
notes the magnification of ray i for sources at redshift 2, 
and lQL(i;2) = 1 if ray i is of type II with lens impact 
parameter ^ 0.5 arcsec, and zero otherwise, 



+ x 



-0 



dx , 



and 



"'^(2) = L|,™(2)/L;(2). The factor |At(i;2) 



(11) 



Eq. ( 10 1 accounts for the fact that we assume a uniform dis- 
tribution of sources in the source plane, while our ray sam- 
pling method uniformly samples the image plane. The ratio 
of the cumulative luminosity distributions quantifies the in- 
crease of source counts due to a lower luminosity threshold 
for detection in regions of higher magnification. 

Under the assumption that the total number of quasars 
is not affected by magnification bias, the number of lensed 



quasars A^ql in the sample of Inada et al. ( 2008 1 is obtained 
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by integrating over the observed redshift distribution nq (z) 
of the quasar sampl^™ 

iVQL 



J weighted by the cross-section: 

roc 

/ nQ{z)aQ-L{z)dz . (12) 

"'0 



If we consider lensing by dark matter alone, we predict 
less than one lensed quasar with image splitting 1 arcsec. 
For lensing by both luminous and dark matter, three lensed 
quasars with image splitting 1 arcsec are predicted. This is 
still small compared to the 11 lensed quasars in the sample. 
One reason for this discrepancy could be that our simple ray- 
selection criterion Iql does not provide a very good approxi- 
mation for the fraction of strongly lensed sources with image 
splittings > 1 arcsec]^ Furthermore, we certainly underes- 
timated the magnification bias[^ Other possible reasons for 
the discrepancy may be effects of non-spherical baryon dis- 
tribution or halo contraction due to baryon condensation. 
These are not considered here. Imitating the contraction ef- 
fect by simply doubling the stellar mass in the galaxies, we 
predict about 8 lensed quasars, much closer to the observed 
numbers of lensed quasars. These results indicate that pre- 
dictions are very sensitive to detailed assumptions about 
galaxy and dark matter structure on small scales, and are 
thus unlikely to provide robust constraints on cosmological 
parameters. 



4 SUMMARY 

In this work, we have studied how the stellar components 
of galaxies affect predictions for gravitational lensing in the 
concordance ACDM cosmology. Our results were obtained 
by shooting random rays though a series of lens planes cre- 
ated from the Millennium Simulation. The dark-matter com- 
ponent on the lens planes was constructed directly from the 
dark-matter particle distribution of the simulation, while the 
distribution of the luminous matter was obtained from semi- 
analytic galaxy formation models run on stored merger trees 
from the simulation. 

In Sec. I3.1l we discussed the influence of stellar mass on 
the statistical distribution of point-source magnifications. 
Although this distribution is almost unchanged for magni- 
fications i-L ~ 1, the galaxies induce a noticeable increase of 
the probability for very high and very low magnifications. 

In Sec. |3.2| we presented optical depths for images of 
small sources that are highly magnified, strongly distorted 
or belong to multiply imaged sources. We compared the re- 
sults obtained by including both dark and luminous matter 



http://www-utap.phys.s. u-tokyo.ac.jp/~sdss/sqls/ 
The criterion Iql effectively excludes (the abundant) lenses 
with rE < 0.5 arcsec, which do not produce image splittings 
> 0.5 arcsec. Undesirably, the criterion also excludes those multi- 
image systems with image separations > 0.5 arcsec where the 
type-II image is near the centre of the lens. 

The magnitude limit for the quasar sample is icor ^ 19.1, 
where icor is the Galactic-extinction corrected magnitude of the 
brightest image with an image-separation dependent admixture 
of the fainter image (s) and the lens galaxy jOguri et al.|[2006| . 
Moreover, the brightest image is often the primary image, which 
is of type I, but not of type II. 



to those obtained for dark matter alone. We find that the in- 
clusion of the luminous matter greatly enhances the strong- 
lensing optical depths compared to a 'dark-matter only' 
universe. Our results in Sec. 13.31 show that the increase is 
partly due to a significant enhancement of the strong-lensing 
cross-section of group and cluster halos with virial masses 
M > 10^^h~^ Mq. In addition, the stellar matter leads to 
significant strong lensing in smaller halos, which would not 
cause noticeable strong lensing otherwise. Although these 
halos have typical Einstein radii < 1 arcsec, their large 
abundance outweighs their small individual cross-sections, 
leading to a bimodal distribution of integrated cross-sections 
with halo mass both for strongly distorted images and for 
multiply imaged sources at redshifts z^ < 1. 



In Sec. |3.4| we estimated optical depths for strongly 
lensed images formed at larger distance from centre of their 
lens. Although the optical depths for lensing by both the 
dark and luminous matter are much larger than for lens- 
ing by dark matter alone for images at radii < 1 arcsec, for 
images at larger radii the optical depths are much more sim- 
ilar. For images at radii ^ 5 arcsec, the optical depths differ 
by at most a factor of two. For radii > 10 arcsec, there is 
almost no enhancement due to the galaxies, in agreement 
with earlier studies (e.g. [Meneghetti et al.|2003| [Dalai et al.| 
|2004) . 

Our results are consistent with the observed radial 
distribution of multiply imaged quasars with splittings ^ 
1 arcsec in the recent SDSS quasar sample of |Inada et al.| 
(2008) only if the effect of the galaxies is taken into ac- 
count. The total number of lensed quasars predicted by our 
standard model is still quite low compared to the number 
observed, indicating a need to include the effects of bary- 
onic dissipation on the dark matter distribution in order to 
explain the data fully. 

An obvious extension to the work presented here will be 
ray-tracing over finite fields to study the effects of lensing 
on sources with finite extent. With a realistic distribution 
of source properties, the results could be compared directly 
with observations of massive galaxy clusters, where many 
distorted images and multi-image systems have been studied 
in some detail ( Broadhurst et al.|2005 Halkola et al.|2006 |. 

Here, we used only the stellar mass and the size in- 
formation from the semi-analytic galaxy formation models. 
In future work, the morphology and luminosity information 
could be added to simulate galaxy-galaxy lensing surveys 
with given selection functions. This will also allow us to take 
into account the effects of the light of lens galaxies on the 
detectability of strongly lensed images near the lens centre 
(Meneghetti et al. 2008) . With some improvements in the 
modelling of the galaxies, in particular by assuming realistic 
elliptical mass and light profiles for the stellar components, 
such simulated surveys will provide detailed predictions for 
galaxy-galaxy lensing. 

The dark-matter distribution in galaxies, groups and 
clusters does not merely provide an arena for physical pro- 
cesses in the baryonic gas, it is also subject to modification 
by these processes ( [Blumenthal et ar]|1986[ [Gnedin et al.| 



2004 1 . Baryon condensation increases the dark-matter den- 



sity in the inner parts of halos, and so affects their dark- 
matter lensing properties ( [Puchwein et al.||2005| |Jing et ah] 
2006 Rozo et al.|2006 1. The semi-analytic models used here 
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do not yet describe these processes, but they should be in- 
cluded in future work aiming at higher precision. 
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